Maackia amurensis seed lectin (MASL) and soluble human podoplanin (shPDPN) sequence analysis and effects on human oral squamous cell carcinoma (OSCC) cell migration and viability.

Maackia amurensis lectins serve as research and botanical agents that bind to sialic residues on proteins. For example, M. amurensis seed lectin (MASL) targets the sialic acid modified podoplanin (PDPN) receptor to suppress arthritic chondrocyte inflammation, and inhibit tumor cell growth and motility. However, M. amurensis lectin nomenclature and composition are not clearly defined. Here, we sought to definitively characterize MASL and its effects on tumor cell behavior. We utilized SDS-PAGE and LC-MS/MS to find that M. amurensis lectins can be divided into two groups. MASL is a member of one group which is composed of subunits that form dimers, evidently mediated by a cysteine residue in the carboxy region of the protein. In contrast to MASL, members of the other group do not dimerize under nonreducing conditions. These data also indicate that MASL is composed of 4 isoforms with an identical amino acid sequence, but unique glycosylation sites. We also produced a novel recombinant soluble human PDPN receptor (shPDPN) with 17 threonine residues glycosylated with sialic acid moieties with potential to act as a ligand trap that inhibits OSCC cell growth and motility. In addition, we report here that MASL targets PDPN with very strong binding kinetics in the nanomolar range. Moreover, we confirm that MASL can inhibit the growth and motility of human oral squamous cell carcinoma (OSCC) cells that express the PDPN receptor. Taken together, these data characterize M. amurensis lectins into two major groups based on their intrinsic properties, clarify the composition of MASL and its subunit isoform sequence and glycosylation sites, define sialic acid modifications on the PDPN receptor and its ability to act as a ligand trap, quantitate MASL binding to PDPN with KD in the nanomolar range, and verify the ability of MASL to serve as a potential anticancer agent.

Benchmarking long-read aligners and SV callers for structural variation detection in Oxford nanopore sequencing data.

Structural variants (SVs) are one of the significant types of DNA mutations and are typically defined as larger-than-50-bp genomic alterations that include insertions, deletions, duplications, inversions, and translocations. These modifications can profoundly impact the phenotypic characteristics and contribute to disorders like cancer, response to treatment, and infections. Four long-read aligners and five SV callers have been evaluated using three Oxford Nanopore NGS human genome datasets in terms of precision, recall, and F1-score statistical metrics, depth of coverage, and speed of analysis. The best SV caller regarding recall, precision, and F1-score when matched with different aligners at different coverage levels tend to vary depending on the dataset and the specific SV types being analyzed. However, based on our findings, Sniffles and CuteSV tend to perform well across different aligners and coverage levels, followed by SVIM, PBSV, and SVDSS in the last place. The CuteSV caller has the highest average F1-score (82.51%) and recall (78.50%), and Sniffles has the highest average precision value (94.33%). Minimap2 as an aligner and Sniffles as an SV caller act as a strong base for the pipeline of SV calling because of their high speed and reasonable accomplishment. PBSV has a lower average F1-score, precision, and recall and may generate more false positives and overlook some actual SVs. Our results are valuable in the comprehensive evaluation of popular SV callers and aligners as they provide insight into the performance of several long-read aligners and SV callers and serve as a reference for researchers in selecting the most suitable tools for SV detection.

Multi-locus sequence analysis unveils a novel genus of filarial nematodes associated with ticks in French Guiana.

Filarial nematodes of the Dipetalonema lineage include tick-borne filarioids that infect both domestic and wild vertebrate hosts, but they remain understudied in many cases. In this study, we conducted a molecular characterization of a Dipetalonema-like filarioid (DLF) recently identified in two tick species in French Guiana, South America. While the cox1 mitochondrial gene was the sole marker initially sequenced for describing DLF, its classification and phylogenetic relationship with other members of the Dipetalonema lineage were unclear. Therefore, we better characterized DLF through the sequencing of six additional gene markers and conducted phylogenetic analyses. Based on this multi-locus typing scheme, DLF exhibited significant divergence from known genera and species of filarioids, or other sequences available in public databases, suggesting its potential classification as a novel genus within the Dipetalonema lineage. Phylogenetic analyses further unveiled a close evolutionary relationship between DLF and all other filarioids associated with Acari (ticks and mites) within a robust monophyletic subclade in the Dipetalonema lineage. Overall, these findings confirm the existence of a specialized, Acari-borne group of filarioids and underscore the need for comprehensive investigations into their epidemiology and potential impact on animal health.

Title: Une analyse de séquences multi-locus dévoile un nouveau genre de nématodes filaires, associé aux tiques en Guyane française. Abstract: Les filaires de la lignée Dipetalonema comprennent des espèces transmises par les tiques qui infectent à la fois des hôtes vertébrés domestiques et sauvages, mais qui restent sous-étudiées dans de nombreux cas. Dans cette étude, nous avons réalisé une caractérisation moléculaire d'un filarioïde ressemblant à Dipetalonema (FRD) récemment identifié dans deux espèces de tiques en Guyane française. Alors que la séquence du gène mitochondrial cox1 était le seul marqueur génétique initialement séquencé pour décrire FRD, sa classification et sa relation phylogénétique avec d'autres membres de la lignée Dipetalonema étaient incertaines. Par conséquent, nous avons caractérisé plus précisément DLF en séquençant six gènes supplémentaires et en réalisant des analyses phylogénétiques. Sur la base de ce typage multi-locus, FRD présentait une divergence significative par rapport aux genres et espèces connus de filarioïdes, ou à d'autres séquences disponibles dans les bases de données publiques, suggérant sa classification potentielle en tant que nouveau genre au sein de la lignée Dipetalonema. Les analyses phylogénétiques ont en outre révélé une relation évolutive étroite entre FRD et tous les autres filarioïdes associés aux tiques et acariens au sein d'un sous-clade monophylétique dans la lignée Dipetalonema. Dans l'ensemble, ces résultats confirment l'existence d'un groupe spécialisé de filarioïdes transmis par les tiques et acariens et soulignent la nécessité d'études approfondies sur leur épidémiologie et leur impact potentiel sur la santé animale.

Genomic characterization of equine influenza A subtype H3N8 viruses by long read sequencing and functional analyses of the PB1-F2 virulence factor of A/equine/Paris/1/2018.

Equine influenza virus (EIV) remains a threat to horses, despite the availability of vaccines. Strategies to monitor the virus and prevent potential vaccine failure revolve around serological assays, RT-qPCR amplification, and sequencing the viral hemagglutinin (HA) and neuraminidase (NA) genes. These approaches overlook the contribution of other viral proteins in driving virulence. This study assesses the potential of long-read nanopore sequencing for fast and precise sequencing of circulating equine influenza viruses. Therefore, two French Florida Clade 1 strains, including the one circulating in winter 2018-2019 exhibiting more pronounced pathogenicity than usual, as well as the two currently OIE-recommended vaccine strains, were sequenced. Our results demonstrated the reliability of this sequencing method in generating accurate sequences. Sequence analysis of HA revealed a subtle antigenic drift in the French EIV strains, with specific substitutions, such as T163I in A/equine/Paris/1/2018 and the N188T mutation in post-2015 strains; both substitutions were in antigenic site B. Antigenic site E exhibited modifications in post-2018 strains, with the N63D substitution. Segment 2 sequencing also revealed that the A/equine/Paris/1/2018 strain encodes a longer variant of the PB1-F2 protein when compared to other Florida clade 1 strains (90 amino acids long versus 81 amino acids long). Further biological and biochemistry assays demonstrated that this PB1-F2 variant has enhanced abilities to abolish the mitochondrial membrane potential ΔΨm and permeabilize synthetic membranes. Altogether, our results highlight the interest in rapidly characterizing the complete genome of circulating strains with next-generation sequencing technologies to adapt vaccines and identify specific virulence markers of EIV.

Full genome sequence analysis of African swine fever virus isolates from Cameroon.

African swine fever (ASF) is a devastating disease of domestic pigs that has spread across the globe since its introduction into Georgia in 2007. The etiological agent is a large double-stranded DNA virus with a genome of 170 to 180 kb in length depending on the isolate. Much of the differences in genome length between isolates are due to variations in the copy number of five different multigene families that are encoded in repetitive regions that are towards the termini of the covalently closed ends of the genome. Molecular epidemiology of African swine fever virus (ASFV) is primarily based on Sanger sequencing of a few conserved and variable regions, but due to the stability of the dsDNA genome changes in the variable regions occur relatively slowly. Observations in Europe and Asia have shown that changes in other genetic loci can occur and that this could be useful in molecular tracking. ASFV has been circulating in Western Africa for at least forty years. It is therefore reasonable to assume that changes may have accumulated in regions of the genome other than the standard targets over the years. At present only one full genome sequence is available for an isolate from Western Africa, that of a highly virulent isolate collected from Benin during an outbreak in 1997. In Cameroon, ASFV was first reported in 1981 and outbreaks have been reported to the present day and is considered endemic. Here we report three full genome sequences from Cameroon isolates of 1982, 1994 and 2018 outbreaks and identify novel single nucleotide polymorphisms and insertion-deletions that may prove useful for molecular epidemiology studies in Western Africa and beyond.

Sequence Analysis of Six Candidate Genes in Miniature Schnauzers with Primary Hypertriglyceridemia.

Miniature Schnauzers are predisposed to primary hypertriglyceridemia (HTG). In this study, we performed whole genome sequencing (WGS) of eight Miniature Schnauzers with primary HTG and screened for risk variants in six HTG candidate genes: LPL, APOC2, APOA5, GPIHBP1, LMF1, and APOE. Variants were filtered to identify those present in ≥2 Miniature Schnauzers with primary HTG and uncommon (<10% allele frequency) in a WGS variant database including 613 dogs from 61 other breeds. Three variants passed filtering: an APOE TATA box deletion, an LMF1 intronic SNP, and a GPIHBP1 missense variant. The APOE and GPIHBP1 variants were genotyped in a cohort of 108 Miniature Schnauzers, including 68 with primary HTG and 40 controls. A multivariable regression model, including age and sex, did not identify an effect of APOE (estimate = 0.18, std. error = 0.14; p = 0.20) or GPIHBP1 genotypes (estimate = -0.26, std. error = 0.42; p = 0.54) on triglyceride concentration. In conclusion, we did not identify a monogenic cause for primary HTG in Miniature Schnauzers in the six genes evaluated. However, if HTG in Miniature Schnauzers is a complex disease resulting from the cumulative effects of multiple variants and environment, the identified variants cannot be ruled out as contributing factors.

Genome Sequence Analysis of Native Xenorhabdus Strains Isolated from Entomopathogenic Nematodes in Argentina.

Entomopathogenic nematodes from the genus Steinernema (Nematoda: Steinernematidae) are capable of causing the rapid killing of insect hosts, facilitated by their association with symbiotic Gram-negative bacteria in the genus Xenorhabdus (Enterobacterales: Morganellaceae), positioning them as interesting candidate tools for the control of insect pests. In spite of this, only a limited number of species from this bacterial genus have been identified from their nematode hosts and their insecticidal properties documented. This study aimed to perform the genome sequence analysis of fourteen Xenorhabdus strains that were isolated from Steinernema nematodes in Argentina. All of the strains were found to be able of killing 7th instar larvae of Galleria mellonella (L.) (Lepidoptera: Pyralidae). Their sequenced genomes harbour 110 putative insecticidal proteins including Tc, Txp, Mcf, Pra/Prb and App homologs, plus other virulence factors such as putative nematocidal proteins, chitinases and secondary metabolite gene clusters for the synthesis of different bioactive compounds. Maximum-likelihood phylogenetic analysis plus average nucleotide identity calculations strongly suggested that three strains should be considered novel species. The species name for strains PSL and Reich (same species according to % ANI) is proposed as Xenorhabdus littoralis sp. nov., whereas strain 12 is proposed as Xenorhabdus santafensis sp. nov. In this work, we present a dual insight into the biocidal potential and diversity of the Xenorhabdus genus, demonstrated by different numbers of putative insecticidal genes and biosynthetic gene clusters, along with a fresh exploration of the species within this genus.

Passionfruit Genomic Database (PGD): a comprehensive resource for passionfruit genomics.

Passionfruit (Passiflora edulis) is a significant fruit crop in the commercial sector, owing to its high nutritional and medicinal value. The advent of high-throughput genomics sequencing technology has led to the publication of a vast amount of passionfruit omics data, encompassing complete genome sequences and transcriptome data under diverse stress conditions. To facilitate the efficient integration, storage, and analysis of these large-scale datasets, and to enable researchers to effectively utilize these omics data, we developed the first passionfruit genome database (PGD). The PGD platform comprises a diverse range of functional modules, including a genome browser, search function, heatmap, gene expression patterns, various tools, sequence alignment, and batch download, thereby providing a user-friendly interface. Additionally, supplementary practical tools have been developed for the PGD, such as gene family analysis tools, gene ontology (GO) terms, a pathway enrichment analysis, and other data analysis and mining tools, which enhance the data's utilization value. By leveraging the database's robust scalability, the intention is to continue to collect and integrate passionfruit omics data in the PGD, providing comprehensive and in-depth support for passionfruit research. The PGD is freely accessible via http://passionfruit.com.cn .

Dynamics of Low-Level Viremia and Immune Activation after Switching to a Darunavir-Based Regimen.

There is an ongoing debate regarding whether low-level viremia (LLV), in particular persistent LLV, during HIV treatment with optimal adherence originates from low-level viral replication, viral production, or both. We performed an observational study in 30 individuals with LLV who switched to a boosted darunavir (DRV)-based therapy. In-depth virological analyses were used to characterize the viral population and the (activity) of the viral reservoir. Immune activation was examined using cell-bound and soluble markers. The primary outcome was defined as the effect on HIV-RNA and was categorized by responders (<50 cp/mL) or non-responders (>50 cp/mL). At week 24, 53% of the individuals were considered responders, 40% non-responders, and 7% could not be assigned. Sequencing showed no evolution or selection of drug resistance in the non-responders. Production of defective virus with mutations in either the protease (D25N) or RT active site contributed to persistent LLV in two individuals. We show that in about half of the study participants, the switch to a DRV-based regimen resulted in a viral response indicative of ongoing low-level viral replication as the cause of LLV before the switch. Our data confirm that in clinical management, high genetic barrier drugs like DRV are a safe choice, irrespective of the source of LLV.

Whole genome-wide sequence analysis of long-lived families (Long-Life Family Study) identifies MTUS2 gene associated with late-onset Alzheimer's disease.

INTRODUCTION: Late-onset Alzheimer's disease (LOAD) has a strong genetic component. Participants in Long-Life Family Study (LLFS) exhibit delayed onset of dementia, offering a unique opportunity to investigate LOAD genetics. METHODS: We conducted a whole genome sequence analysis of 3475 LLFS members. Genetic associations were examined in six independent studies (N = 14,260) with a wide range of LOAD risk. Association analysis in a sub-sample of the LLFS cohort (N = 1739) evaluated the association of LOAD variants with beta amyloid (Aß) levels. RESULTS: We identified several single nucleotide polymorphisms (SNPs) in tight linkage disequilibrium within the MTUS2 gene associated with LOAD (rs73154407, p = 7.6 × 10-9). Association of MTUS2 variants with LOAD was observed in the five independent studies and was significantly stronger within high levels of Aß42/40 ratio compared to lower amyloid. DISCUSSION: MTUS2 encodes a microtubule associated protein implicated in the development and function of the nervous system, making it a plausible candidate to investigate LOAD biology. HIGHLIGHTS: Long-Life Family Study (LLFS) families may harbor late onset Alzheimer's dementia (LOAD) variants. LLFS whole genome sequence analysis identified MTUS2 gene variants associated with LOAD. The observed LLFS variants generalized to cohorts with wide range of LOAD risk. The association of MTUS2 with LOAD was stronger within high levels of beta amyloid. Our results provide evidence for MTUS2 gene as a novel LOAD candidate locus.

Complete genome sequence analysis of Bacillus velezensis A5, a promising biocontrol agent from the Pacific Ocean.

Tobacco bacterial wilt (TBW) caused by Ralstonia solanacearum is a serious soil-borne disease, which seriously damages the growth of tobacco crops. Bacillus velezensis A5 was isolated from 3000 m deep-sea sediments of the Pacific Ocean, and was found to be antagonistic to TBW. Here, we report the complete genome sequence of strain A5, which has a 4,000,699-bp single circular chromosome with 3827 genes and a G + C content of 46.44%, 87 tRNAs, and 27 rRNAs. A total of 12 gene clusters were identified in the genome of strain A5, which were responsible for the biosynthesis of antibacterial compounds, including surfactin, bacillaene, fengycin, difficidin, bacillibactin, and bacilysin. Additionally, strain A5 was found to contain a series of genes related to the biosynthesis of carbohydrate-active enzymes and secreted proteins. Our results indicate that strain A5 can be considered a promising biocontrol agent against TBW in agricultural fields.

Green Waste Compost Impacts Microbial Functions Related to Carbohydrate Use and Active Dispersal in Plant Pathogen-Infested Soil.

The effects of compost on physical and chemical characteristics of soil are well-studied but impacts on soil microbiomes are poorly understood. This research tested effects of green waste compost on bacterial communities in soil infested with the plant pathogen Fusarium oxysporum. Compost was added to pathogen-infested soil and maintained in mesocosms in a greenhouse experiment and replicated growth chamber experiments. Bacteria and F. oxysporum abundance were quantified using quantitative PCR. Taxonomic and functional characteristics of bacterial communities were measured using shotgun metagenome sequencing. Compost significantly increased bacterial abundance 8 weeks after amendment in one experiment. Compost increased concentrations of chemical characteristics of soil, including phosphorus, potassium, organic matter, and pH. In all experiments, compost significantly reduced abundance of F. oxysporum and altered the taxonomic composition of soil bacterial communities. Sixteen bacterial genera were significantly increased from compost in every experiment, potentially playing a role in pathogen suppression. In all experiments, there was a consistent negative effect of compost on functions related to carbohydrate use and a positive effect on bacteria with flagella. Results from this work demonstrate that compost can reduce the abundance of soilborne plant pathogens and raise questions about the role of microbes in plant pathogen suppression.

ClusTrast: a short read de novo transcript isoform assembler guided by clustered contigs.

BACKGROUND: Transcriptome assembly from RNA-sequencing data in species without a reliable reference genome has to be performed de novo, but studies have shown that de novo methods often have inadequate ability to reconstruct transcript isoforms. We address this issue by constructing an assembly pipeline whose main purpose is to produce a comprehensive set of transcript isoforms. RESULTS: We present the de novo transcript isoform assembler ClusTrast, which takes short read RNA-seq data as input, assembles a primary assembly, clusters a set of guiding contigs, aligns the short reads to the guiding contigs, assembles each clustered set of short reads individually, and merges the primary and clusterwise assemblies into the final assembly. We tested ClusTrast on real datasets from six eukaryotic species, and showed that ClusTrast reconstructed more expressed known isoforms than any of the other tested de novo assemblers, at a moderate reduction in precision. For recall, ClusTrast was on top in the lower end of expression levels (<15% percentile) for all tested datasets, and over the entire range for almost all datasets. Reference transcripts were often (35-69% for the six datasets) reconstructed to at least 95% of their length by ClusTrast, and more than half of reference transcripts (58-81%) were reconstructed with contigs that exhibited polymorphism, measuring on a subset of reliably predicted contigs. ClusTrast recall increased when using a union of assembled transcripts from more than one assembly tool as primary assembly. CONCLUSION: We suggest that ClusTrast can be a useful tool for studying isoforms in species without a reliable reference genome, in particular when the goal is to produce a comprehensive transcriptome set with polymorphic variants.

Sequence analysis of the Spike, RNA-dependent RNA polymerase, and protease genes reveals a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia.

During the Covid-19 pandemic, the resurgence of SARS-CoV-2 was due to the development of novel variants of concern (VOC). Thus, genomic surveillance is essential to monitor continuing evolution of SARS-CoV-2 and to track the emergence of novel variants. In this study, we performed phylogenetic, mutation, and selection pressure analyses of the Spike, nsp12, nsp3, and nsp5 genes of SARS-CoV-2 isolates circulating in Yogyakarta and Central Java provinces, Indonesia from May 2021 to February 2022. Various bioinformatics tools were employed to investigate the evolutionary dynamics of distinct SARS-CoV-2 isolates. During the study period, 213 and 139 isolates of Omicron and Delta variants were identified, respectively. Particularly in the Spike gene, mutations were significantly more abundant in Omicron than in Delta variants. Consistently, in all of four genes studied, the substitution rates of Omicron were higher than that of Delta variants, especially in the Spike and nsp12 genes. In addition, selective pressure analysis revealed several sites that were positively selected in particular genes, implying that these sites were functionally essential for virus evolution. In conclusion, our study demonstrated a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing.

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

The Site/Group Extended Data Format and Tools.

Comparative sequence analysis permits unraveling the molecular processes underlying gene evolution. Many statistical methods generate candidate positions within genes, such as fast or slowly evolving sites, coevolving groups of residues, sites undergoing positive selection, or changes in evolutionary rates. Understanding the functional causes of these evolutionary patterns requires combining the results of these analyses and mapping them onto molecular structures, a complex task involving distinct coordinate referential systems. To ease this task, we introduce the site/group extended data format, a simple text format to store (groups of) site annotations. We developed a toolset, the SgedTools, which permits site/group extended data file manipulation, creating them from various software outputs and translating coordinates between individual sequences, alignments, and three-dimensional structures. The package also includes a Monte-Carlo procedure to generate random site samples, possibly conditioning on site-specific features. This eases the statistical testing of evolutionary hypotheses, accounting for the structural properties of the encoded molecules.

Multi-locus sequence analysis reveals phylogenetically segregated Entamoeba histolytica population.

Amoebiasis, caused by the enteric parasite, Entamoeba histolytica, is one of the major food- and water-borne parasitic diseases in developing countries with improper sanitation and poor hygiene. Infection with E. histolytica has diverse disease outcomes, which are determined by the genetic diversity of the infecting strains. Comparative genetic analysis of infecting E. histolytica strains associated with differential disease outcomes from different geographical regions of the world is important to identify the specific genetic patterns of the pathogen that trigger certain disease outcomes of Amoebiasis. The strategy is able to elucidate the genealogical relation and population structure of infecting E. histolytica strains from different geographical regions. In the present study, we have performed a comparative genetic analysis of circulating E. histolytica strains identified from different parts of the world, including our study region, based on five tRNA-linked short tandem repeat (STR) loci (i.e., D-A, NK2, R-R, STGA-D and A-L) and evaluated their potential associations with differential disease outcomes of Amoebiasis. A number of regional-specific, emerging haplotypes of E. histolytica, significantly associated with specific disease outcomes have been identified. Haplotypes, which have a significant positive association with asymptomatic and amoebic liver abscess outcomes, showed a significant negative association with diarrheal outcome, or vice versa. Comparative multi-locus analysis revealed that E. histolytica isolates from our study region are phylogenetically segregated from the isolates of other geographical regions. This study provides a crucial overview of the population structure and emerging pattern of the enteric parasite, E. histolytica.

Small polymorphisms are a source of ancestral bias in structural variant breakpoint placement.

High-quality genome assemblies and sophisticated algorithms have increased sensitivity for a wide range of variant types, and breakpoint accuracy for structural variants (SVs, ≥50 bp) has improved to near base pair precision. Despite these advances, many SV breakpoint locations are subject to systematic bias affecting variant representation. To understand why SV breakpoints are inconsistent across samples, we reanalyzed 64 phased haplotypes constructed from long-read assemblies released by the Human Genome Structural Variation Consortium (HGSVC). We identify 882 SV insertions and 180 SV deletions with variable breakpoints not anchored in tandem repeats (TRs) or segmental duplications (SDs). SVs called from aligned sequencing reads increase breakpoint disagreements by 2×-16×. Sequence accuracy had a minimal impact on breakpoints, but we observe a strong effect of ancestry. We confirm that SNP and indel polymorphisms are enriched at shifted breakpoints and are also absent from variant callsets. Breakpoint homology increases the likelihood of imprecise SV calls and the distance they are shifted, and tandem duplications are the most heavily affected SVs. Because graph genome methods normalize SV calls across samples, we investigated graphs generated by two different methods and find the resulting breakpoints are subject to other technical biases affecting breakpoint accuracy. The breakpoint inconsistencies we characterize affect ∼5% of the SVs called in a human genome and can impact variant interpretation and annotation. These limitations underscore a need for algorithm development to improve SV databases, mitigate the impact of ancestry on breakpoints, and increase the value of callsets for investigating breakpoint features.

Using molecular methods to delineate norovirus outbreaks: a systematic review.

Noroviruses are among the major causative agents of human acute gastroenteritis, and the nature of norovirus outbreaks can differ considerably. The number of single-nucleotide polymorphisms (SNPs) between strains is used to assess their relationships. There is currently no universally accepted cutoff value for clustering strains that define an outbreak or linking the individuals involved. This study was conducted to estimate the threshold value of genomic variations among related strains within norovirus outbreaks. We carried out a literature search in the PubMed and Web of Science databases. SNP rates were defined as the number of SNPs/sequence length (bp) × 100%. The Mann-Whitney U-test was used in comparisons of the distribution of SNP rates for different sequence regions, genogroups (GI and GII), transmission routes, and sequencing methods. A total of 25 articles reporting on 108 norovirus outbreaks were included. In 99.1% of the outbreaks, the SNP rates were below 0.50%, and in 89.8%, the SNP rates were under 0.20%. Outbreak strains showed higher SNP rates when the P2 domain was used for sequence analysis (Z = -2.652, p = 0.008) and when an NGS method was used (Z = -3.686, p < 0.001). Outbreaks caused by different norovirus genotypes showed no significant difference in SNP rates. Compared with person-to-person outbreaks, SNP rates were lower in common-source outbreaks, but no significant difference was found when differences in sequencing methods were taken into consideraton. SNP rates under 0.20% and 0.50% could be considered as the rigorous and relaxed threshold, respectively, of strain similarity within a norovirus outbreak. More data are needed to evaluate differences within and between various norovirus outbreaks.

Common Features of Environmental Mycobacterium chelonae from Colorado Using Partial and Whole Genomic Sequence Analyses.

Nontuberculous mycobacteria (NTM) are environmentally acquired opportunistic pathogens that cause chronic lung disease in susceptible individuals. While presumed to be ubiquitous in built and natural environments, NTM environmental studies are limited. While environmental sampling campaigns have been performed in geographic areas of high NTM disease burden, NTM species diversity is less defined among areas of lower disease burden like Colorado. In Colorado, metals such as molybdenum have been correlated with increased risk for NTM infection, yet environmental NTM species diversity has not yet been widely studied. Based on prior regression modeling, three areas of predicted high, moderate, and low NTM risk were identified for environmental sampling in Colorado. Ice, plumbing biofilms, and sink tap water samples were collected from publicly accessible freshwater sources. All samples were microbiologically cultured and NTM were identified using partial rpoB gene sequencing. From these samples, areas of moderate risk were more likely to be NTM positive. NTM recovery from ice was more common than recovery from plumbing biofilms or tap water. Overall, nine different NTM species were identified, including clinically important Mycobacterium chelonae. MinION technology was used to whole genome sequence and compare mutational differences between six M. chelonae genomes, representing three environmental isolates from this study and three other M. chelonae isolates from other sources. Drug resistance genes and prophages were common findings among environmentally derived M. chelonae, promoting the need for expanded environmental sampling campaigns to improve our current understanding of NTM species abundance while opening new avenues for improved targeted drug therapies.